Analysis of a tuneable coupler for superconducting phase qubits 
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This paper presents a theoretical analysis of the recently realized tuneable coupler for supercon- 
ducting phase qubits (R. C. Bialczak et al., Ref . [T^ . The coupling can be turned off by compensating 
a negative mutual inductance with a tuneable Josephson inductance. The main coupling in this sys- 
tem is of the XX type and can be zeroed exactly, while there is also a small undesired contribution 
of the ZZ type. We calculate both couplings as functions of the tuning parameter (bias current) 
and focus on the residual coupling in the OFF regime. In particular, we show that for typical 
experimental parameters the coupling OFF/ON ratio is few times 10"'^, and it may be zeroed by 
proper choice of parameters. The remaining errors due to physical presence of the coupler are on 
the order of 10"^. 

PACS numbers: 03.67.Lx, 85.25. Cp 



I. INTRODUCTION 

Superconducting qubits^ are potential building blocks 
of a quantum computer. Among their advantages in com- 
parison with other qubit types are an efficient control 
with voltage/current/microwave pulses and use of a well- 
developed technology suitable for large scale integration. 
Recent demonstrations of simple quantum algorithms''"^ 
and three-qubit entanglement "^"^ with superconducting 
qubits are important steps towards a practical quantum 
computation. 

In the standard idea of a gate-based quantum 
computation^ it is important that the qubits are decou- 
pled from each other for most of the time, and the cou- 
pling of a desired type between two (or three) qubits is 
switched on only for a quantum gate operation, while it is 
switched off again after that. Since the superconducting 
qubits cannot be physically moved in space, such cou- 
pling/decoupling should be realized by changing control 
parameters of a circuit. The simplest idea is to tune the 
qubits in resonance with each other for efficient coupling 
and move them out of resonance for decoupling (see., e.g., 
Refs.Q-H)- However, this requires avoiding unwanted res- 
onances, and with increasing number of qubits may lead 
to the problem of "spectral crowding" . Even more im- 
portant limitation of this approach is the following. Be- 
cause the effective coupling strength when the qubits are 
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detuned by an energy Ai? is of order j t^E, where g 
is the tuned value, the ratio of the "switched off' and 
"switched on" coupling strengths is OFF/ON ~ .g/Ai?. 
This "OFF/ON ratio" characterizes a coupler's ability to 
successfully turn on and off the coupling between qubits. 
l\E/h is limited to about a GHz in current supercon- 
ducting qubit devices. To make the ratio small then re- 
quires g to be small, which makes gate operations slow. 
For example, to realize an OFF/ON ratio of 10^^ when 
/S.E/h — 1 GHz would require g to be IMHz, which is 
unacceptably small. 

A different idea is to introduce an extra element be- 
tween the qubits: an adjustable coupler, which can turn 
the coupling on and off. It is a much better approach 
from the architecture point of view since it allows easier 
design of complex quantum circuits. This has motivated 
several experimental^"^' and theoretical^'^"''^ studies of 
adjustable couplers for superconducting qubits. 

In this paper we theoretically analyze operation of a 
recently realized^- tuneable coupler for superconduct- 
ing phase qubits, which demonstrated current-controlled 
tuning of the XX type coupling from to 100 MHz. 
In the next Section we discuss the Hamiltonian of the 
analyzed Josephson circuit and our definition of the two- 
qubit coupling frequencies ^xx and ^zz , corresponding 
to the XX and ZZ types of interaction (for our system 
typically |rixx| ^ \^zz\, so the main coupling is of the 
XX type). In Sec. HI we find flxx and fizz in a sim- 
ple semiclassical way, while in Sec. IV similar results are 
obtained in the lowest-order quantum analysis. We show 
that both rixx and ^zz can cross zero as functions of the 
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FIG. 1: The analyzed scheme of two-qubit coupling, which is 
controlled by the bias current Ib of the coupling Josephson 
junction. Indices 1 and 2 refer to the two qubits, and index 
3 refers to the coupling junction. Coupling inductors L4 and 
L5 have a negative mutual inductance —M. The current Ib 
controls the Josephson inductance L3 of the coupling junction, 
which effectively adds to — M. An additional small coupling 
is via capacitance Ca- 



control parameter (bias current), but typically not simul- 
taneously, thus leading to a non-vanishing residual cou- 
pling, which is discussed in Sec. V. In particular, we show 
that typical OFF/ ON ratio for the coupler is few times 
10"*^; however, a minor modification of the experimental 
circuilji^ (addition of a small coupling capacitance) can 
zero the residual coupling, thus zeroing the OFF/ON ra- 
tio. Actually, this does not mean complete decoupling, 
because we rely on the two-qubit description of a more 
complicated circuit. The remaining coupling effects are 
also discussed in Sec. V and are shown to lead to errors 
on the order of 10~^. Section VI presents numerical re- 
sults of the quantum analysis; for typical experimental 
parameters they are close to the analytical results. Sec- 
tion VII is the Conclusion. In Appendix we discuss the 
position and momentum matrix elements for an oscilla- 
tor with weak cubic nonlinearity, and derive improved 
analytics for flxx and 0,zz- 



II. SYSTEM AND HAMILTONIAN 

Let us consider the system^'^^ shown in Fig. [Tl which 
consists of two flux-biased phase qubits^l characterized 
by capacitances Ci and C2, inductances Li and L2, and 
Josephson energies Eji and Ej2 of the junctions (in 
Fig. [1] the superconducting phases across these Josephson 
junctions are denoted as (pi and (^2)- The qubits are cou- 
pled via an additional Josephson junction (characterized 
by C3 and Ej^) with an adjustable bias current Ib', the 
qubits are connected to this junction via inductances L4 
and L5, which have a negative mutual inductance —M 
(M > 0). We also introduce the qubit coupling via a 
very small capacitance Ca, which was not implemented 
in the experiment but may be important in future ex- 
periments for turning the coupling off more precisely. 

The general idea of this schemei^ is that the Joseph- 
son inductance L3 of the middle junction is essentially in 
series with the mutual inductance — M, and therefore (in 



absence of Ca) the coupling is expected to be (crudely) 
proportional to L3 — M. Then varying L3 by varying the 
bias current Is, it is possible to adjust the qubit-qubit 
coupling strength, which is expected to cross zero when 
L3 « M. 

The Hamiltonian of the system can be derived in the 
standard wa y^^i^^ and written in terms of the phases 
(i = 1,2,3) across the three Josephson junctions and 
the conjugated momenta {[4>i,pj] = ihSij), which are 
the corresponding node charges, multiplied by $0 = 
$o/27r = h/2e: 
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where Ci = Ci + ^2^/(^2 + Ca) and C2 = C2 + 
CiCa/{Ci + Ca) are the effective qubit capacitances, 
C3 = C3 (introduced for notational convenience), 0e^i 
and (f)e,2 are the external dimensionless qubit fluxes, 0e,4 
and 0e 5 are the external dimensionless fluxes through 
the loops containing L4 and L5 (we will often assume 

0e,4 = (^e,5 = 0), Ca = {l/Ca + 1/Ci + 1/C2)-^ is the ef- 
fective coupling capacitance, and renormalized coupling 
inductances are 



L4 _ L5 _ M _ ^ M2 
L4 L5 M LjiL^ ' 



(2) 



All terms in Eq. ([T]) have clear physical meaning*22 Intro- 
ducing the shifted variables = (pi — 4'i,st, where the 
set {4>i^st\ corresponds to the minimum of the potential 
energy, we rewrite the Hamiltonian ^ as 



H = Hi+H2+H3 + Hi 



int ; 
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z = 1,2,3, 
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where the potentials Ui have minima at Scpi — 0, and the 
corresponding plasma frequencies w^^pi = {LiCi)"^^'^ are 
governed by the effective inductances 
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where in the last equation we expressed Ej^ cos 03. st in 
terms of the corresponding junction current I^ st — Ib 
and the critical current I^^ci = Ej^/^q. 

The Hamiltonians Hi and H2 correspond to the sep- 
arated qubits; in absence of coupling two lowest eigen- 
states in each of them correspond to the logic states 
|0) and Notice that because of an anharmonicity 
of the potentials, the qubit frequencies ioi and uj2 (de- 
fined via energy difference between |1) and |0) for un- 
coupled qubits) are slightly smaller than the plasma fre- 
quencies uJi.pl and 0^2, pi- The coupler is characterized 
by the Hamiltonian H^; its similarly defined frequency 
wa is slightly smaller than wa^pi. In the experiment^ 0-13 
was almost an order of magnitude higher than u}i^2, and 
the coupler had only virtual excitations. In our analysis 
we also assume absence of real excitations in the cou- 
pler; however, in general we do not assume W3 ^ aji.2 
(except specially mentioned), we only assume absence of 
resonance between the coupler and the qubits. 

Our goal is to calculate the coupling between the two 
qubits due to Hmt (we assume a weak coupling) . In gen- 
eral, the coupling between two logic qubits can be char- 
acterized by 9 real parameters. However, since for phase 
qubits the energies of the states |0) and |1) are signif- 
icantly different, we can use the rotating wave approx- 
imation (RWA), i.e. neglect terms creating two excita- 
tions or annihilating two excitations; then there are only 
3 coupling parameters in the rotating frame: 

+ ^ (.«4^) - a^y%^^) 4- (9) 

where superscripts of the Pauli matrices indicate qubit 
numbering. Moreover, from the symmetry of the inter- 
action Hamiltonian ([5]) it follows that ftxY = (because 
matrix elements of 6(j)i are real, and for pi they are imagi- 
nary) ; therefore our goal is to calculate only two coupling 
frequencies: ftxx and fizz- Notice that since the XX 
and YY interactions are indistinguishable in the RWA 
[both correspond to a^_:^^ a^_^^ + a^_^^a^_^^], we may rewrite 

the first term in Eq. ^ as {hflxx/2) cr^V^\ 

The considered system of Fig. [T] has a large Hilbert 
space, which can in principle be reduced to a two-qubit 
space in a variety of ways, giving in general different val- 
ues of flxx and fizz- To avoid ambiguity, we define the 
coupling frequencies in terms of the exact eigenstates of 
the full physical system. In particular, we associate the 
two-qubit logic states |00) and |11) with the correspond- 
ing eigenstates of the full Hamiltomian ([3]) and denote 
their energies as Eqq and En] for this association we 
start with the product-state using the ground state for 
the coupler, and then find the nearest eigenstate. Simi- 
larly, instead of trying to define the uncoupled logic states 
|01) and 1 10) (that is ambiguous, though a natural defini- 
tion can be based on the dressed states discussed in Sec. 
IV) , we deal with the eigenstates resulting from their cou- 



pling, which are associated with two exact eigenstates of 
the full system; their energies are denoted as E+ and 
{E+ > E^). 

Then the coupling fizz is defined as 

hflzz = Eoo + En-iE+ + E^), (10) 

which is obviously consistent with Eq. (|9]) for logic qubits 
(a similar definition has been used in Rcf. 2). The cou- 
pling flxx can be defined as the minimal splitting in the 
avoided level crossing between |01) and |10), i.e. as 

\hflxx\ = min {E+-E^), (11) 

with the sign of flxx easily obtained by comparing with 
Eq. ([9]). Actually, the definition (ITlT) cannot be applied 
to an arbitrary qubit detuning wi — W2 (in a symmetric 
case it works only for degenerate qubits); such gener- 
alization of flxx definition can be done by comparing 
exact eigenstates with the standard avoided level cross- 
ing behavior (discussed in more detail in Sees. IV and V). 
Notice that our definitions of flxx and fizz do not need 
any assumption of a weak coupling (this is their main 
advantage); however, a weak coupling will be assumed in 
derivation of analytical results. 

III. SIMPLE SEMICLASSICAL ANALYSIS 

Let us first calculate flxx and fizz in a simple, es- 
sentially electrical engineering way (we will see later that 
the result is close to the quantum result). For simplic- 
ity in this section we assume (L4, L5) >• (Zi, L2, £3, ill), 
Ca <C (Ci,C2), so that the coupling is weak and the 
tilde signs in many cases can be avoided. We also re- 
place the middle junction with the effective inductance 

L3 = ^ol^llcr-llst ~ ^o/^llcr-ll this ap- 
proximation L3 w L3). 

The coupling flxx corresponds to the frequency split- 
ting between the symmetric and antisymmetric modes of 
the two-qubit oscillations. So, let us assume degenerate 
qubits, wi = W2 = Wqb, and find the splitting in the clas- 
sical linear system. Notice that at frequency Wqb the ca- 
pacitance C3 is equivalent to the inductance — I/^lo'^^C^), 
and therefore the parallel connection of L3 and C3 is 
equivalent to the inductance 

3 l~u;l^L,Cs l-(^qbM)2' ^ ' 

notice that here wa — a;3.pi for the coupler since we as- 
sume a linear system. 

Suppose u! « cjqb is a classical eigenfrequency, and 
the first qubit voltage is Vie*"*. Then using the phasor 
representation, we find the current through L4 as I4 = 
Vil {iojL^); it induces the voltage V^p = iw(L|'^ — M)/4 in 
the coupling inductances, which causes the current = 
T4p/(iwL5) = Vi{Lf - M)/{iLLiLiLz) flowing through 
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is into the second qubit. Adding this current to the cur- 
rent la = iujCaVi through C^, we get the total current 
I2 = Ib+Ia, flowing into the second qubit. The extra cur- 
rent I2 is equivalent to changing the qubit capacitance C2 
by AC2 = -l2/{iujV2), where V2 = ±Viy^CJC^ is the 
second qubit voltage for the symmetric and antisymmet- 
ric modes (the factor -y/Ci/Z^ comes from the condition 
of equal energies in the two qubits). The effective change 
of the capacitance slightly changes the oscillation fre- 
quency {L2C2)~^^^ , so the eigenfrequency can be found 
as w = Wqb(l — AC2/2C2). Therefore, the frequency 
splitting due to coupling is jilxxl = i^qb|AC2|/C2, and 
substituting AC2 we finally find 



IV. QUANTUM ANALYSIS (ANALYTICS) 

For the quantum analysis let us rewrite the interaction 
Hamiltonian ([SJ as 

Hint = -^^13(01 + al){a3 + 4) + K23{a2 + 4)(«3 + 4) 
+Ki2{ai + a\){a2 + 4) + ^72(01 - 4)(«2 - 4): (15) 

where Ui + aj — S4>iy^2miUJi/h, — a| = ipi^2/hmiUJi, 
and 

Ki3 = ~ fci3, K23 = 5 K23, (16) 



Af-L3/[l-(u;qbM)'] , CaOJ. 



qb 



L4_L5Wqb-\/ C1C2 



V C1C2 



(13) 



M „ Ca 

Ki2^^^ki2, i^ia = -777^77— : 

L4L5 C162 4fci2 



kij — 



where the explicit expression (IT2|) for has been used, 
and the sign of ^xx is determined by noticing that a 
positive ftxx should make the frequency (energy) of the 
symmetric mode larger than for the antisymmetric mode 
[see Eq. (HI)]. 

The most important observation is that flxx depends 
on the bias current Ib, which changes L3, and for a 
proper biasing the coupling flxx can be zeroed exactly. 
If the correction due to the C^-term is small and also 
<C 1 (as in the experiments^), then ftxx is ze- 
roed when L3 w M. 

The coupling hQzzcrz^ a'"^^ /A in Eq. © originates 
from anharmonicity of the qubit potentials and cor- 
responding difference between the average Josephson 
phases for states |1) and |0), which we denote as A0^q 



2y/miUJimjUjj 



(17) 
(18) 



(we use the creation/annihilation operators a| and 
only for brevity of notatons; in their normalization we 
use the frequency uJi between two lowest eigenstates in- 
stead of the plasma frequency) . 

In order to find flxx , we have to solve the Schrodinger 
equation H\ip±) = E±\i/j±) for the two eingenstates |-0±) 
{E+ > E-), corresponding to the coupled logic states 
1 10) and |01). The wavefunction can be written in the 
product-state basis as \tlj±) = a±|100) -I- /3±, |001) -I- 
where in this notation we show the energy levels \'nin3n2) 
of the first qubit, the coupling oscillator (in the middle), 
and the second qubit, and the terms not shown explic- 
itly should be relatively small in the weak coupling case. 
Comparing the amplitudes and energies with the stan- 
dard avoided level crossing behavior, let us define two 



for the ith qubit. This leads to the extra dc current coupling frequencies, V.^^ and ftxx^ 
$o^0ioV-^4 through the inductance L4, when the first 



^10 

qubit changes state from |0) and |1), and similar dc cur- 
rent change ^oA(I)\q /L^ through L5 for the second qubit. 
As a result, the state |11) acquires an additional magnetic 
interaction energy, which is the product of these two cur- 
rents multiplied^! by M — L3. This corresponds to 



- 



E+ — Ezi 



1 + {a±/P±? 



(19) 



zz 



(14) 



We have to define two frequencies because in general 
«+//?+ ^ —j3-/a-, in contrast to the ideal case of two 
logic qubits; this is the price to pay when the two-qubit 
language is applied to a more complicated physical sys- 
tem. However, the difference l^^^x ~ ^xjc I typically 
very small; moreover, for degenerate qubits in a symmet- 



ric system 



-I- _ 

XX 



^xx (because then |q;±//3±| — 1), 



Comparing this result with Eq. (|13l) for flxx, we see 
absence of the contribution due to Ca (which is small 
anyway) and a similar proportionality to M — L3, though 
without the correction 1 — (a;qb/w3)^. This means that 
by changing the bias current Ib (which affects L3), the 
coupling fizz can be zeroed, and this happens close to 
the point where flxx is zeroed. Except for the vicinity 
of the crossing point, \il,zz/^xx\ ^ 1 because A0^q is 
small (compared to the ground state width) for a weak 
anharmonicity; therefore flxx is the main coupling in 
our system. 



so that we need only a single frequency flxx, which in 
this case coincides with the definition (ITT|) . We will ne- 
glect the difference between fi^x ^^d ^Ixx unless spe- 
cially mentioned (the difference is important in the case 
of strongly detuned qubits). 

For the analysis it is convenient to express a solu- 
tion of the Schrodinger equation Hlip) = E\^) as \4>) = 
'^IV'ioo) + /^IV'ooi); where the dressed states |'0ioo) ^-^d 
IV'ooi) ^'"S defined in the following way. The state IV'ioo) 
expanded in the product-state basis has the contribu- 
tion from the state |100) with amplitude 1 and zero 
contribution from the state |001), i.e., (V'loollOO) ~ 1 
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and (■0ioo|OOl) = 0. Also, IV'ioo) satisfies equation 
(n|-ff IV'ioo) = -^('^IV'ioo) fo^' basis elements \n) = 
\n1_nzn2) except |100) and |001). The dressed state 
IV'ooi) is defined similarly, except now (tAqq-^IIOO) = 
and (■0001 1001) — 1- Notice that a dressed state is not 
a solution of an eigenvalue problem; for a given energy 

it is a solution of an inhomogeneous systems of linear 
equations. Also notice that we do not need to normalize 
the wavefunctions. 

Constructing the dressed states in this way, we have to 
satisfy (self-consistently for E) only two remaining equa- 
tions to solve the Schrodinger equation: 

a(100|i?|V?5o> + /?(100|iJ|Ci> = Ea, (20) 
a(001|i/|Vf5o> + /?(001|i?K^i) = Ep. (21) 

Using linear algebra it is easy to prove the reciprocity 
relation (100|i?|V'o(5i) = (001|ff |V'i5o)* (in our case the 
complex conjugation is actually not needed since the 
matrix elements are real), and therefore Eqs. and 
pT|) are similar to the standard equations for an avoided 



level crossing. Hence, the matrix elements Ef^^ = 
(100|77|i/;5'go) and E^^^ = {mi\H\i)%^) play the role of 
renormalized self-energies of the two-qubit logic states 
1 10) and |01), while the two-qubit coupling can be calcu- 
lated as 

^xx = ^(OOl|i?i„t|0?5o) = ^(100|ffintK\), (22) 

so that the eigenenergies are given by the usual formula 

^^n^]/2 [in 



E± = [E% + E^^, ± J(£;d5o - E^^.y 



Eq. (|22p we wrote Hint instead of H because there is 
obviously no contribution from the non-interacting part] . 
Notice that the dressed states and therefore the matrix 
elements depend on energy E, in contrast to the standard 
level crossing. This leads to a slight difference of ilxx for 
the eigenstates i?+ and E-, and also makes calculations 
using Eq. slightly different from the definition p^ . 

To find ^xx analytically, let us assume that the three 
oscillators are linear. Hi = hjji{a\ai 4- 1/2), and use the 
lowest-order perturbation theory. In the first order 



l^fSo) = |100)+i^i3 



|010) ^/2|210) 



E — epio E — £210 



,K.,J^,^K.,,Kt,,^. ,23) 



E~eii 



E — €201 



where the energies e of the basis states are only due to 
non-interacting part Hi + H2 + H3 of the Hamiltonian 
dS]). Then from Eq. ([22]) we obtain 



nxx - 1{Ki2-K^2)+Iki3K23 (-^ + 

(24) 

For degenerate qubits and weak coupling we can use ap- 
proximation E £100 — eooi in this equation, so that 
E - eoio « -^(^^3 - Wqb) and E - em sa -h{uj3 + Wqb); 



then using explicit expressions ([16 
elements we finally obtain 



^xx = 



A/ - L*/[l - (c^qbM)'] , C>qb 



L^L^Uab \/ C1C2 



TSl) for the matrix 



(25) 



V C1C2 



where L*^ = i3(l + f )(1 + |f ), CI = ^((71^2/^1^2), 
and we used 0^3 — [^36*3]^^/^. Comparing this equation 
with the classical result ([13]), we see that the results coin- 
cide under assumptions used for the classical derivation. 

Validity of the perturbation theory requires assump- 
tions \Ki'i/h{LOz± UJcih)\ < 1, |iC23/^(w3 ± Wqb)| < 1, 

|ifi2/Sa;qb| ^ 1, and jiffj/^qbl ^ 1, which basically 
mean that in Eq. (|25l) the contributions to 0.xx due to 
M and C* should be much smaller than Wqb, and the 

contribution L\lL^L^LOa^\J C1C2 should be much smaller 
than |a;3 — ti;qbp/w3. Notice that we do not need an as- 
sumption ajqb/w3 <C 1, we only need absence of resonance 
between these frequencies; in fact, W3 can be even smaller 
than Wqb. 

To find analytics for ^zz, it is necessary to consider 
nonlinear oscillators, because this is what we expect from 
the quasiclassical analysis and also because in the lin- 
ear case for degenerate qubits eioi — £200 = £002, and 
therefore there is an ambiguity in defining the logic state 
1 11). In order to calculate ^zz via Eq. ([TU)) . we need to 
find the eigenenergies i?oo and En, while calculation of 



E+^E. 



Efnn + Enm has been already discussed above. 



100 "T -^001 

To find the ground state energy Eqo we introduce the 
dressed state IV'ooo) as a state satisfying the Schrodinger 
equation {n\H\ip) = E{n\'ilj) (for a given E) for all basis 
elements |n) except |000), and also satisfying condition 
(OOOIV'too) = 1- In a similar way as above, we construct 
iV'ooo) in the first order perturbation theory and then 
find _Eoo as i?ooo = (000|i/ I'f/'ooo)' using approximation 
E ~ £000 in the construction of the dressed state. To 
find the eigenenergy En, we introduce the dressed state 
IV'ioi) in a similar way, then calculate it in the first or- 
der, and then find i?ii as Efgi = (101|77|i/)ioi) assuming 
E « £101 for the dressed state. 

Even though this is a straightforward procedure, now 
there are infinitely many terms in the first-order dressed 
states because of the nonlinearity, and there are still 
many terms even if we keep only lowest orders in non- 
linearity. However, most of the contributions to the 
energies cancel each other in the combination Mlzz ~ 
Eqqq + E^Qi — Ef^Q — Eq^^ , and the largest non-canceling 
contributions yield 



nzz = 



6162 



K12 



13^^23 



£101 — Cm 



where 6, is defined for zth oscillator as 



(1|50,|1)-(O|50,|O) 



(26) 



(27) 



Notice that fizz depends on the nonlinearity of qubits 
(via bi and 62), while nonlinearity of the coupling junc- 
tion gives only a small correction (see Appendix) to the 
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second term in Eq. (|26p , which is neglected in the lowest 
order. Also notice that in Eq. (^5)) we neglected terms 
proportional to Kf2, because they are on the same order 
as the neglected terms ~ {Ki3K23)^ ■ Using the defini- 
tions (|16p - p^ . we rewrite ((26)) in the form 



A I 



zz 



(28) 



2 L4Lr-,^UJlU}2\/CiC; 

which coincides with the classical result under as- 
sumptions used for the classical result, since A(/)^q = 

In deriving Eq. p6| we have used the lowest order 
of the perturbation theory. However, there are higher- 
order terms, which are significantly enhanced because 
the basis state |101) is close to resonance with the states 
|200) and |002) even for degenerate qubits. Let us ac- 
count for this effect by analyzing the repulsion between 
these levels and computing the corresponding shift of 
the eigenenergy En. Following the above formalism for 
^xx, we find the level splitting due to interaction be- 
tween |101) and |200) to be 5'|ii),|20) = 2(101|7?int|V'20o)- 
Then writing IV'200) same way as in Eq, 



/,dr ^ 

we find >5'|ii)j20) ~ V^fi^xx, which is essentially the 
same result as for a qubit interacting with a resonator 3^ 
Because of the level repulsion, the eigenenergy En has 
a shift by [e200 - eioi ± \/ (^200 - eioi)^ + '^fi^^xx)/'^' 
which in the dispersive case cioi — £200 ^ l^^jfjc | becomes 
^^^xx/2(eioi - £200) (we assume eioi > £200, so the 
shift is up in energy). In another notation eioi — £200 — 
h{ui2 — ijJi+Suii), where by hSuJi we denote the correction 
for the second excited level energy, 2ei — eo — £2, for ith 
qubit. A similar shift up in energy for En comes from 
the interaction with the level |002). Adding these two 
contributions, we modify Eq. ((28|) to become 
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zz 



(72 



1 



1 



UJ2 — i^l + SuJi U>i — UJ2 + SuJ2 



(29) 



Notice that for fizz we in general consider different qubit 
frequencies uJi and uj2, while in Eq. (1^ for ^xx we 
assumed nearly degenerate qubits; however, unless qubit 
detuning significantly affects ilxx (that will be discussed 
later), we can use definition ujqh = ('^i+'^2)/2 in Eq. p5|). 

Since Qzz has a major dependence on the qubit non- 
linearity, let us discuss it in more detail (see also Ap- 
pendix). For ith oscillator potential with an additional 
cubic term, it is convenient to characterize nonlinearity 
by the ratio Ni = J/bar.i/^i^i.pi, where J7bar,i is the barrier 
height (assumed to be at 6(j)j > 0), so that Ni is crudely 
the number of levels in the quantum well (iVi^2 ~ 5 in 
typical experiments with phase qubits^). For a weak cu- 
bic nonlinearity {Ni ^ 1) one can derive^ the following 
approximations: 



(5/36iV,) Lj, 



(30) 



Therefore, away from the point where flxx ~ 0, and 
neglecting corrections due to non-zero C^, Wqb/ws, and 
W2 — , the ratio of couplings is 



zz 



1 



l8iNi+N2)n 



XX 



n 



XX 



6VN^ 



(31) 



which is quite small for typical experimental parameters. 
Notice that for Ni — 5 (which is typically used for qubits) 
the numerical values bi = 0.289 and SoJi/ui — 0.0378 
are significantly different from what is expected from the 
large-iV analytics pop (in the cubic approximation for 
the qubit potential bi and SoJi/uJi depend only on A^i). 



V. RESIDUAL COUPLING 

Both ^xx and Qzz may cross zero when L3 is var- 
ied by adjusting the bias current Ib- However, they are 
typically zeroed at different values of L3, that prevents 
turning the two-qubit coupling completely off. Since 
\^zz/^xx\ ^ 1 away from the zero-crossing points, let 
us characterize the residual coupling by the Z Z-cowpling 
value ri'^l at the point where ^xx = 0. Using Eqs. (^5)) 
and (|29l) . we find for degenerate qubits 



zz 



bib2 
2 



C: 



2 t^qb 



\fChC2 



(32) 



where we assumed Wqb/w3 ^ 1, as in the experiment. 

It is natural to characterize the OFF/ON ratio for the 
adjustable coupling by the ratio 57^|/f2^^, where f^xx 
is the "fully on" AA-coupling. Let us assume the op- 
erating regime in which the ON-coupling corresponds to 
zero bias current;^ and at this point M / ~ 2 [see Eq. 

m] , 

an estimate 



Then l^ON 



M /2L4L5Ujqi,\/ C1C2, and we find 



(-)rcs 
"ZZ 



/n 



ON 
XX 



1"2 



'■'qb 



7^3) 



(33) 



In particular, for Ca — 0, Ni — N2 — 5, and 
~ 1/5 (typical experimental parameters), this 
gives OFF/ON ~ 3 x lO'^. The capacitance Ca ^ 
M / L^L^Lol needed to zero the OFF/ON ratio is then 
around 0.6 fF for typical experimental parameters L4 5 ~ 
3 nH, M ~ 200 pH, 0J3/2n ~ 30 GHz. 

Actually, since our analytics is only the leading-order 
calculation, while in fi^z have an almost exact can- 
celation of contributions due to M and L3, we cannot 
expect that Eqs. ([32]) and ([33]) are accurate even in the 
leading order. Nevertheless, we can trust the result that 
the OFF / ON ratio is typically quite small, because both 
6162 ~ NiN2)~^ and the terms in the brackets in 
Eq. ([551) are small. Moreover, the OFF/ON ratio can be 
made exactly zero by choosing proper values for Wqb/ws 
and Ca- Even if our analytics missed a small term in 
Eq. ([551) . the OFF/ON ratio can be zeroed either by in- 
creasing Ca or decreasing 0^3, since this moves iV^zz 
the opposite directions. 
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Since f2^^/i7^^ is very small or even zero, we have 
to carefully consider other effects, which do not vanish 
when both ftxx and fizz discussed above are zero. One 
of such effects becomes clear when we consider the case 
of a strong qubit detuning, \u!i — u!2\ ^ \^xx\- In this 
case one would expect that the eigenstate close to the 
state 1 100) should have a negligible contribution of the 
state |001) and vise versa (so that the logic states |10) 
and 1 01) are decoupled); however, actually these con- 
tributions cannot be decreased to zero. As seen from 
Eq. (|24p. ^xx depends on energy E, and therefore it 
is slightly different for the two eigenstates with energies 
E+ and The difference Aflxx = ^xx ~ ^xx 

approximately —4:Ki3K23{E+ — E^)/h?u)1, which can be 
rewritten as 



Af7 



XX 



\UJI - i02\ 



-2n 



ON I'^l ~ ^2] 



XX' 



1C2 W3 

where we assumed uj^ 3> 2 and used 



(34) 

.3 ^ M for 

nearly OFF coupling (while M/L^ ~ 2 in the ON regime 
at Ib =0). This means that if we zero the amplitude of 
the state |100) in one of the eigenstates, there will still be 
a non-zero amplitude of the state |001) in the other eigen- 
state (in contrast to an ideal two-qubit situation). Choos- 
ing the smallest flxx coupling as ±Ar2xx/2, we obtain 
contributions — i7^^/2cj3 of the wrong states in both 
eigenstates. So, the error occupation is ~ (r2°^/2w3)^, 
which for typical parameters is around 10^®. 

Another effect, which is related to inaccuracy of the 
RWA approximation, can be characterized by the contri- 
bution of the state |101) in the ground state. Using the 
second-order perturbation theory for the dressed state 
IV'ooo) we find 



K12 + Kf^ 



K13K2 



£000 — eioi 



(35) 



which is approximately —flxx/4:'^qb when Qxx is 
not close to zero, exactly as expected for the non- 
RWA contribution {h^xx /'2.)<7+^ cr^+^ from the term 

(?ii7xx/2)cr^''cr^^ in the two-qubit Hamiltonian. How- 
ever, we are mostly interested in the case flxx = 0; 
then Eq. ([55)1 becomes approximately — f2^^/2aj3, and 
the corresponding error occupation is ~ /2uj3)^ ~ 

10~^, same as for the strong-detuning effect. 



VI. NUMERICAL RESULTS 

To calculate flxx and Qzz numerically, we first find 
the phases {0i,st, </'2,st, 03,st}, which correspond to the 
minimum potential energy in Eq. ([1}, and then find 
the eigenfunctions and eigenenergies of the Hamiltonian 
(l3])-(l5]) using the product-state basis of energy levels in 
the three anharmonic oscillators (in this calculation we 
use the cubic approximation for the oscillator poten- 
tials). After that the coupling fizz is calculated from 
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FIG. 2; The coupling flxx/2TT and flzz/^Tr (multiplied by 
5 for clarity) as functions of the bias current Ib- Solid 
lines show numerical results, while dashed lines show ana- 
lytics using Eqs. (|25[) and (|29l) . The circuit parameters are: 
/i,cr = h.cT = 1.5 ^A, J3,cr = 3 /iA, Ci = C2 = 1 pF, 
C3 = 0.1 pF, Ca = 0, Li = L2 = 0.7nH, L4 = is = 3nH, 
M = 0.2 uH; 0e,4 = <^.,5 = 0, iVi = iV2 = 5 (a;i/27r = 
aj2/27r = 6.59 GHz). The panel (b) is a blow-up of the panel 
(a) near the crossing points at positive Ib- We see that 
iixx = at Ib/Is.ct = 0.759, and the residual coupling 
(square symbol) is nf|/27r = -172 kHz. 



To)) , while for ilxx we cal- 



using the definition 



the eigenenergies using Eq, 
culate two values 

(fT9|) . However, for the case of degenerate qubits, which 
we mostly consider below, there is no difference between 
i^j^j^ and f^j^j^ . 

Figure [2ja) shows ^Ixx and il.zz as functions of 
the bias current Ib for the system with the following 
parameters!^ /i cr = ha- = 1-5 fi A, I3 cr — 3/iA, 
Ci = C2 = IpF, C3 = b.lpF, Ca = 0, ii ^ L2 = 
0.7 nH, L4 ^ L5 ^ 3nH, M = 0.2 nH. We assume 
0e,4 = 0e,5 = for the coupler loops, while the qubit ex- 
ternal fluxes (j)e,i ~ 0e,2 are chosen so that for the qubits 
N\ — N2 — "o; this corresponds to the qubit frequencies 
loi/2t: = W2/27r = 6.59 GHz, which are kept constant 
with changing Ib by the compensating change of exter- 
nal fluxes 4>e,i and (/'e.2- The results for ^Izz are multi- 
plied by 5 for clarity (to become visually comparable to 
^xx)- The solid lines in Fig.[2l[a) show numerical results, 
while dashed lines show the analytics using Eqs. ([25|l and 
One can see that overall the analytics gives a pretty 
good approximation. Smaller numerical value for U,xx 
than in analytics can be partially explained by the cor- 
rections shown in Eqs. (|A.7p and (jA.Sp in Appendix. We 
have checked numerically that the beating frequency of 
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the classical small-amplitude oscillations is close to the 
analytical quantum result for flxx shown in Fig. [2{a), 
with a typical difference on the order of 1 MHz. 

The lines in Fig. [2{a) are not symmetric about Ib = 
because of the current through the coupling junction 
coming from the qubits, which adds to Ib (the curves 
are symmetric about Ib = ~0.122/3.cr; this asymme- 
try could be removed if the qubits are biased with op- 
posite fluxes, 4)e.2 ~ ~4>e,ii SO that the currents from 
the qubits compensate each other). At zero bias (ON- 
coupling) the coupling ilxx/'^^T^ is 34.3 MHz, which (an- 
alytically) comes from 85.8 MHz coupling due to the mu- 
tual inductance —M and compensating —51.1 MHz from 
the inductance L3 of the coupling junction (analytical 
total slightly differs from the numerical result). 

At both positive and negative Ib the coupling Vlxx 
crosses zero because of increase of L3, while ^zz barely 
crosses zero because of the similar increase of L3 and al- 
ways positive contribution from the level repulsion effect 
[see Eq. (l29t ]. Figure [2jb) is a blow-up of [2ja) near 
the crossing points at positive Ib- One can see that 
numerically calculated residual coupling (at VLxx = 0) 
is f^f|/27r = -172 kHz, so that the OFF/ON ratio is 
5 X 10"'^. While this value of f^^^ is on the same or- 
der as expected from the analytics (see dashed lines), it 
has the opposite (negative) sign. This apparently hap- 
pens because corrections to analytics in this case have 
a stronger effect than the effect of the term {oJqh/^^zY 
in Eq. (I?5|) : at the point where VI xx = we have 
ujz/^TT — 38.6 GHz, so (wqb/a;3)^ = 0.029 is really quite 
small (a;3/27r = 49.6 GHz at Ib = 0). Notice that a 
negative value of fi^| makes impossible to zero 51^^ by 
adding the capacitive coupling via Ca- 
in order to make fi^^ positive, we can decrease W3 by 
increasing the coupling junction capacitance C3. Figure 
[3l^a) shows results for the same circuit with C3 — 0.3 pF, 
which is 3 times larger than in Fig. [51 We show only 
vicinity of the crossing points, while the overall shape of 
the curves is quite close to what is shown in Fig. [2Ja); in 
particular, 17xx/27r = 32.5 MHz at Ib = 0. As we see 
from Fig. [3Ja), now 51^1 becomes positive, r2™|/27r = 
49 kHz (at this crossing point ^3/271 ~ 22.9 GHz, and it 
is 28.7 GHz at Ib = 0). The corresponding OFF/ON 
ratio is now 1.5 x 10""^. 

Obviously, ^^zz — ^ for some intermediate value of C3. 
We have calculated that it happens for C3 = 0.253 pF 
(keeping other parameters unchanged). If C3 is larger 
than this value, so that il^z is positive, we can zero the 
residual coupling by adding small coupling capacitance 
Ca- Figure EJb) shows the results for C3 = 0.3 pF and 
Ca = 0.155 fF (other parameters unchanged), in which 
case 17^1 = 0. 

At the point where ^xx — ^zz = 0, we have to pay 
a special attention to other effects which couple the two 
qubits. In particular, we should consider what happens 
when the qubits are detuned in frequency. Figure U shows 
^xx^ ^XX' ^^'^ ^zz (multiplied by 10) as functions of 
the detuning wi — UJ2 for the circuit with C3 = 0.3 pF 




FIG. 3: (a) Same as in Fig. [Sl^b), but for a circuit with 
C3 = 0.3 pF. This makes positive Jlg'|/27r = 49 kHz (square 
symbol), (b) Same as in (a), but with added coupling ca- 
pacitance Ca ~ 0.155 pF, that produces fl^zz = 0, i.e. the 
couplings Qxx and Qzz are zeroed simultaneously. 



and Ca = 0.155 fF (other parameters as above) at the 
bias current Ib — 0.742 /s^ci-; these parameters corre- 
spond to the point ^xx — ^zz = in Fig. ^h). For 
the detuning we change the external fluxes (/)e,i and 0e,2 
(and correspondingly change A''! and N2), while keep- 
ing the frequency ujqh — ('^i + ^^2)72 unchanged (and 
we still assume (peA = 4'e,5 = 0). As seen in Fig. SI 
the couplings ^xx ^^'^ ^xx coincide when uji — uj2] 
however, their difference grows with the qubit detuning. 
The dashed lines show the analytics for zLAflxx /'^ using 
the first expression in Eq. p4|) . There is a significant 
difference between the analytical and numerical results 
because the ratio Wqb/a;3 = 6.59 GHz/23.0 GHz is not 
very small; the next order correction to the analytics by 
the factor [1 -I- 3(a;qb/a;3)^], which comes from Eq. (pi]) . 
accounts for most of the difference. Even though 
grows with the detuning, the corresponding error state 
occupation |a//3|^^ is constant and is only 1.5 x 10~^. 
The detuning also changes ilzz] however, the effect is mi- 
nor, and rizz/'^TT — —42.6 kHz at the detuning of 1 GHz. 
An almost vertical feature on the fizz line at the de- 
tuning of 0.23 GHz is due to the level crossing between 
states 1 101) and |002). It is relatively small because we 
have chosen the operating point with ilxx — in ab- 
sence of detuning; otherwise at the level crossing point 
the coupling V,zz changes by approximately Olxx / 
[see discussion above Eq. d^ ]. 

For the numerical results in this section we have used 
the cubic approximation for the qubit and coupling os- 
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FIG. 4: Solid lines: numerical results for ^xx' ^xx' ^'i'^ 
fizz (multiplied by 10) as functions of the qubit detuning 
uji — UJ2 for the parameters of Fig. |3jb) at the point where 
rixx ~ ^zz ~ 0. Dashed lines show the analytics for 
^xx ~ ±AS7js:x/2 from Eq. (I34|) . Error state occupation 
due to nonzero is 1.5 x 10~^. 



cillator potentials in calculation of the matrix elements 
of the Hamiltonian. We have also done calculations us- 
ing the exact potential and checked that the values of 
flxx and fizz change only slightly, though the residual 
coupling changes more significantly as expected for an 
almost exact cancelation of contributions with opposite 
signs. 



VII. CONCLUSION 

The main goal of this paper has been calculation of the 
two-qubit coupling frequencies ftxx and ilzz for the cir- 
cuit of Fig. [U and analysis of their dependence on the bias 
current Ib , which can be used to turn the coupling on and 
ofFfi^ We have shown that a simple "electrical engineer- 
ing" analysis for i^xx as the beating frequency of two 
classical oscillators gives a result, Eq. which is close 
to the analytical (lowest-order) quantum result, Eq. (pS)) 
(the formulas coincide, except minor renormalization of 
parameters). In turn, the quantum analytics for flxx 
is close to the results of the quantum numerical anal- 
ysis (Fig. [2]); a minor difference is due to higher orders 
in perturbation. The "electrical engineering" analysis for 
^zz, Eq. dH]), is not fully classical; it needs the language 
of quantum energy levels and shows that il,zz originates 
due to anharmonicity of the qubit oscillators, which shifts 
the average flux. However, this analysis corresponds to 
only one term in the quantum analytics, Eq. (1^^ . while 
the other significant contribution is due to the level re- 
pulsion between the states with the single-excitation in 
each qubit and with the double-excitation in one of them. 
The results of the numerical quantum analysis for fizz 
are similar to the quantum analytics (Fig. [2]). 

As expected, our analysis shows that flxx is the main 
two-qubit coupling in the considered circuit, and fizz is 
typically much smaller. Nevertheless, in the analyzed nu- 
merical example using realistic experimental parameters, 
the ratio ftxx/flzz is only around 5 for the coupling 



turned on (small Ib), which means that corrections for 
non-zero fizz in experimental algorithms are necessary. 

The most practically important case is when the cou- 
pling is almost off. The fact that flxx can be zeroed 
exactly is rather trivial: a real number changing sign 
should necessarily cross zero, and flxx obviously changes 
sign when the effect of the coupling junction inductance 
overcompensated the effect of the magnetic coupling 
—M. The coupling ilzz does not necessarily change sign 
because of always positive contribution from the level re- 
pulsion, but fizz should be small when flxx = 0, as 
follows from the analytical result (p9)) . This fact is quite 
beneficial for the ability to turn the coupling almost off, 
and we have defined the residual coupling fl'^^z 
value of fizz when flxx = 0. A natural measure of 
the coupling OFF/ON ratio is fl'^^z/^xx^ where flf^x 
corresponds to the case of small (zero) Ib- Notice that 
this ratio depends on the definition ^ since we compare 
different couplings (for example, if the ZZ-teim was de- 
fined in dni) without the factor 4, then the OFF/ON ratio 
would decrease 4 times). As we found from the analyti- 
cal and numerical calculations, the OFF/ON ratio is few 
times 10"'^ for typical experimental parameters. Most 
importantly, the OFF/ ON ratio can be zeroed exactly 
by properly choosing capacitances C3 and/or Ca- 

Even when the above defined OFF/ON ratio is exactly 
zero, this does not mean that the qubits can be made 
completely decoupled. The reason is that the qubits are 
still physically coupled to the coupling circuit, and the 
discussion of such effects should go beyond the language 
of the coupling of logical qubits, which is characterized 
by only flxx and fizz- In particular, detuning of the 
formally decoupled qubits leads to an erroneous state oc- 
cupation of around (r2^^/2a;3)^, which is on the order 
of 10~^ for typical experimental parameters. Non-RWA 
corrections bring errors of the same order. 

In this paper we sometimes assumed a high fre- 
quency of the coupling oscillator, uj^ ^ Wqb, as in the 
experiment. -'^^ However, the analytical results and 
(l29t . which do not rely on this assumption, show that 
it is not really needed for the operation of the scheme. 
Moreover, the adjustable coupling can be even realized 
for UJ3 < Wqb ; however in that case we have to use positive 
mutual inductance, and we should not expect typically 
small fizz when flxx — 0. 

The calculation of the coupling flxx and fizz in this 
paper is based on the analysis of the eigenfunctions and 
eigenenergies of the whole system, thus avoiding ambi- 
guity of reducing the whole system to two logical qubits. 
However, we have also done the analytical calculations 
by using such reduction. Assuming Wqb/ws ^ 1, we have 
eliminated the coupling junction degree of freedom by ap- 
plying the Schrieffer-Wolf transformation and then pro- 
jecting the resulting Hamiltonian on the coupling junc- 
tion ground state; after that the Hamiltonian has been 
truncated to the two-qubit subspace. The obtained re- 
sults for flxx and fizz basically coincide with Eqs. (pS)) 
and (j28p under the assumptions used. 
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The analyzed tuneable coupler (without Ca) has been 
realized by Bialczak et al.^ and the dependence of the 
coupling frequency ftxx on the bias current Ib has been 
measured experimentally (the coupling ^Izz has not been 
measured). Due to a relatively small critical current /a^cr 
of the coupling junction in the experiment, the coupling 
ilxx is crossing zero at small Ib (see Fig. 4d of Ref. Il2h . 
Therefore, in contrast to the case shown in our Figs. 2 
and 3, the experimental coupler is nearly OFF at Ib — 0. 
Notice that Fig. 4d of Ref. 112| shows —^Ixx (in our nota- 
tion), so it increases with \Ib\, and the dependence on Ib 
is symmetrized by a horizontal shift. Using the experi- 
mental parameters, we have checked that the theoretical 
result for Qxx is quite close to the experimental result 
(Fig. 4d of Ref. shows a fitting by the simple theory, 
which is close to the full theory result). 

For realization of multi-qubit algorithms it is very im- 
portant that the residual coupling ri^z '^^^ zeroed 
by proper design of C3 and Ca- If this is not done, the 
OFF/ON ratio is small (few times 10"'^), but may still be 
significant for complicated algorithms. There is a modifi- 
cation of the scheme of Fig. [1] which may further reduce 
the OFF/ON ratio without using Ca and without pre- 
cise choice of C3. The idea is to add blocking capacitors 
between the qubits and inductors £4,5. Then there will 
be no dc current from the qubits going into the coupling 
circuit, and this will eliminate the classical interaction 
effect leading to ilzz in Eq. (IT4)) . Correspondingly, this 
should eliminate the first term in the quantum result (j29p 
for fizz, so that 57^1 should be very small by itself. In 
the quantum language, this happens because in the mod- 
ified scheme the capacitive interaction between five oscil- 
lators is of the momentum-momentum type (besides one 
phase- phase interaction due to M), and the average mo- 
mentum for any eigenstate of an oscillator is exactly zero. 
From experimental point of view, the scheme with block- 
ing capacitors is more convenient because it eliminates 
the need to adjust external fluxes in the qubits when 
the bias current Ib is changed. We have performed pre- 
liminary quantum calculations for flxx, which confirm 
that flxx crosses zero when the effective inductance 
compensates the magnetic coupling —M, similar to the 
case without blocking capacitors. However, mathemat- 
ically this involves compensation of three dozen quan- 
tum terms of the same order, so we may expect that the 
scheme with blocking capacitors is less robust against 
decoherence than the scheme of Fig. [TJ Such a compar- 
ative analysis of the schemes with and without blocking 
capacitors is a subject of further study. 

This work was supported by NSA and lARPA under 
ARO grant W911NF-08-1-0336. 



Appendix: Corrections due to nonlinearity 

In this appendix we discuss an oscillator with a weak 
cubic nonlinearity, and show next-order corrections for 
^Ixx and fizz due to nonlinearity. 



Let us consider an oscillator with a cubic nonlinearity, 
H = p'^/2m + (mw^j/2)(50)2 - A(50)^, where cjpi is the 
plasma frequency and A > 0, so that there is a finite 
barrier height f/bar — m^tOp^/bAX^ at positive 6(1). It is 
convenient to characterize nonlinearity by the ratio N ~ 
Uhar/f'^pi, SO that N is crudely the number of levels in 
the quantum well. 

Nonlinearity changes the eigenstates eigenenergies 
efc — {k\II\k), and the normalized matrix elements of the 
coordinate and momentum operators, 

(kMl) , , {k\p\l) 



Ckl = Clk 



y^h/2mujpi 



dki — —dik 



(A.l) 

For a weak nonlinearity {N ^ 1) one can derive (similar 
to Ref. [29I ) the following approximations: 

(A.2) 
,(A.3) 




dkk = 



(A.4) 
(A.5) 



N ' V54]V 

where we keep only the lowest-order non-vanishing cor- 
rections, and non-integer numbers are numerical results. 
In particular, these approximations lead to Eq. (1301) for 
b = (cii — coo)\/w/a;pi and 6uj — (2ei — eo ~ £2)/^, where 
a; = (ei — eo)/fi. With next-order corrections, Eq. ([50)) 
becomes 

1 0.28 Suj 5 0.18 

(A.6) 



w 36iV ' 

Now let us discuss corrections for the analytics for flxx 
and flzzi which are next order in nonlinearity, while we 
still use the lowest order in the perturbation theory. This 
modifies Eq. (1^^ to become 



XX 



23 



^01,3 



■-01,3 



^02.3 



■-02.3 



^ — eoio ^ — eiii ^ — £020 -c^ — £121 
where the additional index i = 1,2,3 in Cki,i and dki,i 
is the oscillator number. Now using Eqs. (|A.2|) - ()A.5|) we 
obtain 
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Equation (j26p with account of next-order corrections in 
nonlinearity becomes 



fizz = 



6162 



K 



2i^i3if23C^l,3 2i^i3/^23 C^a 3 



12 



£101 — £111 



£101 — £121 



(A.9) 
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that modifies Eq. [and the first term in Eq. ((291)] to 
become 

f... ^ i±f±S £^^iii±HZM. (A.10) 

A useful technical trick in deriving these results [as well 



as Eq. ([26]) ] is to slightly shift the coordinates S(j)i for the 
three oscillators (by slightly changing {0i,st, '/>2,st, </>3,st}) 
to produce (0|(5(^i|0) — for each oscillator; this signifi- 
cantly reduces the number of terms due to nonlinearity 
in the derivation. 
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